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Cellular senescence is thought to play a major role in age-related diseases, which cause nearly 67% 
of all human deaths worldwide. Recent research in mice showed that exercising mice had higher 
levels of telomerase, an enzyme that helps maintain telomere length, than non-exercising mice. A 
commonly used model for biological aging was proposed by Penna. I propose two modifications of the 
Penna model that incorporate senescence and find analytical steady state solutions following Coe, 
Mao and Cates. I find that models corresponding to delayed senescence have younger populations 
that live longer. 



I. INTRODUCTION 

From early alchemists looking for the elixir of life, to 
modern day researchers, humans have always wanted to 
understand aging. As people age, their cells go through 
replication cycles. Each replication reduces the length 
of the telomeres in the cells. If a cell's telomeres are 
too short, it may not be able to replicate [T]. Cells that 
can no longer replicate are termed senescent. Cellular 
senescence is thought to play a major role in age-related 
diseases, which cause nearly 67% of all human deaths 
worldwide [2] ■ Recent research in mice showed that ex- 
ercising mice had higher levels of telomerase, an enzyme 
that helps maintain telomere length, than non-exercising 
mice[3]. In humans, runners had longer telomeres than 
non-runners [3]. Since shortened telomeres are thought 
to be related to death, this research would seems to in- 
dicate that people who exercise live longer lives. Popula- 
tion studies do indeed show this, however, most studies 
show that only the mean lifespan increases in exercising 
populations, not maximum lifespan [H [5]. 

A commonly used model for biological aging was pro- 
posed by Penna in 19950 [7J- The model looks at death 
from a mutation accumulation standpoint. In the last 10 
years, papers have been published describing methods for 
finding the age distributions created by the Penna model 
without actually simulating the model [8HlO|. However, 
the research on senescence indicates that age-related 
death is not caused by mutation accumulation, but rather 
by an inability to reproduce after a number of cycles. I 
propose two modifications of the Penna model that use 
this mechanism instead of a mutation accumulation, and 
show that a mean lifespan can increase without affecting 
the maximum lifespan, which cannot be done by chang- 
ing parameters in the original Penna model. 

II. CALCULATIONS 

The Penna model assigns a bit-string to each individual 
in the population. Each bit corresponds to a timestep of 
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FIG. 1. The possible paths until death, after reaching I, in 
the SP. 

the simulation. A 1 in the string represents a mutation, 
and a means no mutation. If an individual has gone 
through T Is, then it dies. Each individual can have 
offspring, with probability b. The child's bit-string is 
derived from the the bit-string of the parent, where each 
has probability m of becoming a 1. The Penna model 
ignores positive mutations, because they are rare. The 
length of the bit-string provides a hard limit for lifespan. 
In this paper, I will work with only the T = 1 case. 

In the first proposed modification, which we will call 
the senescent Penna (SP) model, each individual can only 
get one disease, which is essentially the beginning of ag- 
ing. After the individual starts to age, it has probability 
p of staying alive to reproduce at each time step. The 
maximum number of replication cycles is M. This is 
the Penna model, except instead of definitely dying, the 
individual has only a probability of dying. People who 
exercise generally have longer telomeres, so once senes- 
cence starts, they have a smaller probability of dying. 
Higher values of p represent exercising populations. 

The first 1 in an individual's bit string is the age, /, at 
which senescence begins. The number of people alive at 
time j, with age x, is rij(x, I, m) = prij^i(x — \, I, m — 1) if 
x > I, where p is the probability of living after senescence, 
and m is the time since the inception of senescence. If 
x < I, then rij(x, I) — rij-i(x — 1, 1), lib is the probability 
of birth, and e _/3 is the probability of an individual going 
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unmutated, then the number of children born in the next 
time step with disease acquisition time I, Uj+x(0, 1, 0) is 
given by 

rij + i(0, 1,0) = be~ fjl V' nj (x, I, m) 

x m 

+ {l-e-f 3 )be-^Y.T. n ^ l '^ m ) 
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Since m is just the maximum of and x — I, 



(i) 



l-l M+l-l 

^2^2nj(x,l,m) = ^2n(x,l) + ^ n(x,l) (2) 

x m a;— x—l 

where M is the maximum number of replication cy- 
cles allowed after disease. Assuming a steady state, 
n,j + i(x,l) — rij(x,l). Therefore, 
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^2<n(x,l,0) = I -71(0,1,0), 
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x=0 



since each x has the same number of people. Since 
n,j(x, I) = prij_i(x — 1, I) for x > I, and since nj+i(x, I) = 
n,j(x,l) in a steady state 
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Defining as 
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and n(i) = n(0,Z), eq. 1 can be simplified to 
= fcr^nfl) " — + (l-e->e-^ £ »(!') (6) 
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Writing the same equation for I + 1, some algebra leads 
to 

This equation leads to some limiting cases, in order to 
maintain a steady state. Neither the numerator, nor the 
denominator should vanish in eq. 7. 
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However, eq. 7 only gives the time of senescence, not the 
lifespan. In the original Penna model, at I the individual 



dies, but here the individual has only a (1— p) probability 
of dying. The number of individuals who die at age t is 
just 

t 

D(t) = p M n{t — M) + C- - p)p t ~ x n{x), (10) 

x=t-M+l 

assuming t > M. The first term ensures that when t = 
I + M, all the remaining people alive with n(Z) die, not 
just a proportion of 1 — p. 

It is also possible to combine the original Penna model 
with the modification described previously, making a 
combination Penna (CP) model. Some people have a 
genetic tendency for heart disease, or stroke [TT]. These 
diseases are not directly affected by senescence. In a com- 
bination of the two models, there would be two classes 
of death, one instantaneous, as in the original model, 
and one probability based, as in the proposed modifica- 
tion. Then each individual is defined in terms of I and 
lp, where I is the time of a probabilistic mutation and 
lp is the time to a Penna mutation, or instant death. 
nj+i(0, I, lp) can come from rij(x, I' , lp) where V > I and 
lp > lp. Note that I cannot be greater than lp. There- 
fore, 

n j+1 (0,l,l P ) = be- al e-^ p [^2 nj (x,l,l P ) 

X 

+ (l- c -«)(l- e -0)£; £ ^(xj'j'p) 

l'>U'p>lp x 

V>1 x 

+ (l-e- a ) ]T J2 nj (x,l,l'p)], 

i' P >ip x 

(11) 

where e a and e ^ are the probabilities of not getting 
a senescence or death mutation, respectively. The death 
distribution is given by probabilistic and instantaneous 
deaths, 



D(t)= P l ~ X n{x,t) 



=t-M 



(12) 
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+ p M n(t-M, l P ) 

As time passes in the Penna model, the mutations 
slowly move downwards. However, there is no evolution- 
ary pressure preventing an individual with I = from 
reproducing. Since there are no positive mutations, my 
ansatz is that in a steady state, I = for all individuals. 
If lp = 0, the individual does not reproduce, so there 
is an evolutionary pressure for lp > 0. If this happens, 
then any terms summing over different values of I vanish, 
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FIG. 2. Higher values of p result in earlier times of senescence 
in the SP. In this picture, M = 5. 



FIG. 3. The death curve in the SP is a shift of the senescence 
curve, with the amount shifted varying on p. Here M = 5 and 
p = .8. 



so 



n 
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Simplifying in the same manner as in the previous model, 
the steady state form is the same as eq. 7, except 
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However, since I = 0, the time until senescence is always 
0, eq. 12 is not needed. The recursive form just gives the 
age of death. 
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III. RESULTS 

Interestingly, in the SP model higher values of p show 
a "younger population" (fig. 2). With lower p values, an 
individual with low I will not have as many opportunities 
to reproduce, since chances are it will die out soon. This 
provides an evolutionary pressure for higher values of I. 
However, if p is high, then individuals with low I can still 
reproduce. 

While senescence begins later for populations with 
lower p, death comes earlier once senescence is reached. 
The average death time for n(l) is I + (1 — p) + 2p(l — 
p) + 3p 2 (l - p) + . . . + (M + l)p M , where the last term 
has no 1 — p factor, since everybody left alive has to die. 
Higher values of p result in later deaths, as would be ex- 
pected. The later deaths are more apparent after l maxi 



FIG. 4. Higher values of p have a higher proportion of their 
population reach larger ages, but that proportion difference 
is countered in the middle of the death distribution. M = 5 



when any living members of the population are undergo- 
ing senescence. Figures 4 and 5 show the differences in 
percentages, which have to add up to 0. The increased 
proportion of older individuals for higher values of p has 
to be balanced by a reduced proportion of the popula- 
tion at lower ages. The jump in the differences of death 
proportions early in figs. 4 and 5 is caused by M. For 
larger values of p, the 1— p term in eq. 10 is small enough 
to make the p l ~ x term negligible. However, once M is 
reached, the final term has no (1 — p) factor, causing the 
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FIG. 5. This shows the same break as fig. 4, but at age 10, FIG. 7. The death distributions for different values of p 
since M — 10 in the SP here. look approximately the same, but there is a slightly higher 

proportion of people alive at later times for higher values of 

P- 




Time until death 



FIG. 6. The differences in the death distribution curves of 
the SP of M = 5 and M = 10, for different values of p. 



jump. 

The SP model does not show an increase in the mean 
for higher values of p, but it does show that the probabil- 
ity of living to a higher age is greater. Even though the 
probability of living to a high age is greater, the maxi- 
mum age for both populations with higher and lower p is 
still l max +M. Higher values of M push the time of senes- 
cence further forward, since there is less evolutionary in- 
centive for an individual to have a higher I. However, 
just like p, higher values of M also afford a longer time 
until death, balancing out the earlier senescence times. 
Since higher values of M allow for a longer life, the pro- 



FIG. 8. This is the death distribution in the CP, for different 
values of p. Notice that even a slight change in p has a large 
change on the death distribution. 



portion of people alive at a higher age is greater for higher 
values of M. The breaks in fig. 5 are caused by the same 
mechanism as the breaks in figs. 3 and 4, except instead 
of the deaths coming at one specific M, they come at the 
two values of M, 5 and 10. 

The CP model shows a different structure than the SP 
model. It shows a quicker increase in age expectancy, 
and a slower decrease. The combination requires a high 
p, because otherwise q degenerates rapidly to j^— . If p 
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is too low, then the CP model will not be stable without 
a Verhulst factor, 



V = 1 
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where N(t) is the number of people alive at time t, N max 
is the maximum allowed number of people, and V is the 
probability that an individual survives a timestep. This 
also explains why even small differences in the value of p 
show large differences in the age distribution curve. 



IV. CONCLUSION 

In this paper, I showed that simple modifications to the 
Penna model allow for shifts in the lifespan distribution 
without changing the maximum lifespan. Higher values 
of p in the SP result in younger populations, but they die 
later. 

Thinking about the original context for this modifica- 
tion, higher values of p can represent exercising popula- 



tions. The maximum age of two populations with differ- 
ent p and the same M is the same, but the exercising 
population has a higher chance of reaching later ages. 
What is interesting, is that this model shows that the 
exercising population will also be younger. The times 
of senescence for high p are lower than low p, yet the 
probability of reaching a high age is greater. 



The Penna model is a tool used to help us under- 
stand population dynamics. My modifications of the 
Penna model take into account senescence of a popula- 
tion, which is a critical part of the aging process, and help 
to explain the changes in lifespan observed in exercising 
and non-exercising populations. By adjusting e~^,p, and 
M, these should be able to fit actual data. Further im- 
provements to the models could take into account recent 
research in autism which suggests that e - ^ is actually a 
function of time [12], and also looking at positive muta- 
tions. 
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